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High-energy neutrinos, arising from decays of mesons that were produced through the cosmic rays collisions with 
air nuclei, form unavoidable background noise in the astrophysical neutrino detection problem. The atmospheric 
neutrino flux above 1 PeV should be supposedly dominated by the contribution of charmed particle decays. 
These (prompt) neutrinos originated from decays of massive and shortlived particles, , D°, D°, D^ 1 , Aj!~ , 
form the most uncertain fraction of the high-energy atmospheric neutrino flux because of poor explored processes 
of the charm production. Besides, an ambiguity in high-energy behavior of pion and especially kaon production 
cross sections for nucleon-nucleus collisions may affect essentially the calculated neutrino flux. There is the 
energy region where above flux uncertainties superimpose. 

A new calculation presented here reveals sizable differences, up to the factor of 1.8 above 1 TeV, in muon 
neutrino flux predictions obtained with usage of known hadronic models, SIBYLL 2.1 and QGSJET-II. The 
atmospheric neutrino flux in the energy range 10 — 10 7 GeV is computed within ID approach to solve nuclear 
cascade equations in the atmosphere, which takes into account non-scaling behavior of the inclusive cross- 
sections for the particle production, the rise of total inelastic hadron-nucleus cross-sections and nonpower-law 
character of the primary cosmic ray spectrum. This approach was recently tested in the atmospheric muon flux 
calculations pQ. The results of the neutrino flux calculations are compared with the Frejus, AMANDA-II and 
IceCube measurement data. 



1. Introduction 



Atmospheric neutrinos (AN) appear in decays of 
mesons (charged pions, kaons etc.) produced through 
collisions of high-energy cosmic rays with air nuclei. 
The AN flux in the wide energy range remains the is- 
sue of the great interest since the low energy AN flux 
is a research matter in the neutrino oscillations stud- 
ies, and the high energy atmospheric neutrino flux is 
now appearing as the unavoidable background for as- 
trophysical neutrino experiments [lOEHSHllEllSJIH]. 
To present day a lot of AN flux calculations are made, 
among which HaifflHaHaHSIElliailZllIBI (see 
also pjJJ [20] for a review of ID and 3D calculations of 
the AN flux), but so far we don't know, how discrep- 
ancy is strong in the conventional neutrino flux re- 
sulted from various hadronic interaction models, how 
much differences are great due to uncertainties in 
primary cosmic ray spectra and composition in the 
"knee" region. 

In this work we present results of new one- 
dimensional calculation of the atmospheric muon neu- 
trino flux in the range 10-10 7 GeV made with use 
of the hadronic models QGSJET-II 03 [2T], SIBYLL 
2.1 [22] as well as the model by Kimel & Mokhov 
(KM) [23] that were tested also in recent atmospheric 
muon flux calculations [TJ [23]. We compute here 
zenith-angle distribution of the conventional neutrinos 
and compare calculated neutrino energy spectra with 
the data of AMANDA-II and IceCube experiments. 



2. The method and input data 

The calculation is performed on the basis of the 
method [25] of solution of the hadronic cascade equa- 
tions in the atmosphere, which takes into account 
non-scaling behavior of inclusive particle production 
cross-sections, the rise of total inelastic hadron-nuclei 
cross-sections, and the non-power law primary spec- 
trum (see also [TJ[TH1[23]). As the primary cosmic ray 
spectra and composition in wide energy range used 
is the model recently proposed by Zatsepin & Sokol- 
skaya (ZS) j2DJ, which fits well the ATIC-2 experiment 
data [27] and supposedly to be valid up to 100 PeV. 
The ZS proton spectrum at E > 10 6 GeV is com- 
patible with KASCADE data [28] as well the helium 
one is within the range of the KASCADE spectrum 
obtained with the usage of QGSJET 01 and SIBYLL 
models. Alternatively in the energy range 1 — 10 6 GeV 
we use the parameterization by Gaisser, Honda, Lipari 
and Stanev (GH) [2D], the version with the high fit to 
the helium data. Note this version is consistent with 
the data of the KASCADE experiment at E Q > 10 6 
GeV that was obtained (through the EAS simulations) 
with the SIBYLL 2.1. To illustrate the distinction of 
the hadron models employed in the computations, it 
is appropriate to compare the spectrum- weighted mo- 
ments (Table m computed for proton-air interactions 
(for 7 = 1.7): 

l 

ZpJEo) = [ —^dx, (1) 
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Table I Spectrum weighted moments z pc (Eo) calculated 
for 7 = 1.7 



Model 


E , GeV 


Zpp 






Z pn~ 


z pK+ 


z pK~ 




10 2 


0.174 


0.088 


0.043 


0.035 


0.0036 


0.0030 


QGSJET 


10 3 


0.198 


0.094 


0.036 


0.029 


0.0036 


0.0028 


11-03 


10 4 


0.205 


0.090 


0.033 


0.028 


0.0034 


0.0027 




10 2 


0.211 


0.059 


0.036 


0.026 


0.0134 


0.0014 


SIBYLL 


10 3 


0.209 


0.045 


0.038 


0.029 


0.0120 


0.0022 


2.1 


10 4 


0.203 


0.043 


0.037 


0.029 


0.0097 


0.0026 




10 2 


0.178 


0.060 


0.044 


0.027 


0.0051 


0.0015 


KM 


10 3 


0.190 


0.060 


0.046 


0.028 


0.0052 


0.0015 




10 4 


0.182 


0.052 


0.046 


0.029 


0.0052 


0.0015 



where x — E c /Eo, c — p, n, 7r ± ,K ± . The values 
in Table [I] display approximate scaling law both in 
SIBYLL 2.1 and KM and little violation of the scaling 
in the QGSJET-II for p and tt ± . 



3. Atmospheric muon neutrino flux 

Along with major sources of the muon neutri- 
nos, 7r^2 and K^i decays, we consider three-particle 
semileptonic decays, K^ 3 , K^ s , the contribution orig- 
inated from decay chains K — > it — > {K% — > ir + ir~ , 
— > ^ir ), as well as small fraction from the muon 
decays. 

One can neglect the 3D effects in calculations of 
the atmospheric muon neutrino flux near vertical at 
energies E > 1 GeV and at E > 5 GeV in case of 
directions close to horizontal (see [T51 IT5] V 

A comparison of (y^ + v^) flux calculations for 
the three hadronic models under study is made 
in Table [TT1 column 1, 2 and 3 presents the 

fl^io^^Mr, 4? GSJET - n) /^ M) and 
0^ibyll)^(qgsjet-ii) correspondingly, calculated at 
6 = 0° and 90° (in brackets) with usage of the GH 
and ZS primary spectrum. One can see that usage 
of QGSJET-II and SIBYLL models leads to apparent 
difference of the muon neutrino flux, as well as in the 
case of SIBYLL as compared to KM (unlike the muon 
flux, where SIBYLL and KM lead to very similar re- 
sults pQ). On the contrary, the QGSJET-II neutrino 
flux is very close to the KM one: up to 100 TeV the dif- 
ference does not exceed 5% for the GH spectrum and 
10% for the ZS one at = 0°. While the muon flux 
discrepancy in the QGSJET-II and KM predictions is 
about 30% at vertical [1] . The origin of differences is 
evident: the kaon production ambiguity. 

Zenith-angle distributions of the conventional neu- 
trinos, <j> v (E, Q)/<j) v AE, 0°), for the energy range 
1 — 10 5 TeV are shown in Fig. [l] Calculations are 
made with QGSJET-II and SiBYLL 2.1 models both 
for GH and ATIC-2 primary spectra and composition. 



Table II Ratio of the Vf, fluxes at 9 — 0° (90°) calculated 
with the SIBYLL 2.1, QGSJET-II, and KM 



E v , GeV 


1 


2 


3 






GH 




10 2 


1.65 (1.22) 


0.97 (0.85) 


1.65 (1.36) 


10 3 


1.71 (1.46) 


0.96 (0.92) 


1.73 (1.50) 


10 4 


1.60 (1.57) 


0.96 (0.96) 


1.58 (1.55) 


10 5 


1.54 (1.49) 


0.99 (0.96) 


1.46 (1.46) 






ZS 




10 2 


1.58 (1.26) 


1.00 (0.91) 


1.58 (1.38) 


10 3 


1.64 (1.39) 


0.95 (0.92) 


1.73 (1.51) 


10 4 


1.55 (1.46) 


0.96 (0.95) 


1.61 (1.54) 


10 5 


1.37 (1.23) 


0.91 (0.83) 


1.51 (1.48) 


10 6 


1.10 (0.95) 


0.61 (0.55) 


1.80 (1.73) 


10 7 


0.89 (0.75) 


0.48 (0.43) 


1.85 (1.74) 



As was expected, a shape of the angle distribution 
visibly depends on the neutrino energy (at E < 100 
TeV) particularly close to horizontal. The effect of 
hadronic models (as well as of the primary spectrum) 
on the angle distribution is weak. 
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Figure 1: Zenith-angle enhancement of the i/ M + V\± flux. 

Figure [2] shows this work calculations of the neu- 
trino flux (lines) in comparison with the result of Barr, 
Gaisser, Lipari, Robbins and Stanev (BGLRS) 15J) 
obtained with use of the TARGET 2.1 (symbols). All 
these computations are performed for the GH pri- 
mary spectra. As one can see the calculations for KM 
and TARGET 2.1 are in close agreement in the range 
10 - 10 4 GeV (near horizontal) or at E v < 400 GeV 
near vertical. 

The calculation of conventional + fluxes at 
different zenith angles is compared with preliminary 
data of IceCube experiment in Fig. [3] Curves (cos 6 = 
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Figure 2: The two independent calculations for the GH 
spectrum [20] . 



t- 1.0 from top to bottom) display calculated fluxes 
made for GH primary spectra and composition with 
usage of QGSJET-II model, blue points with error 
bars present the IceCube muon neutrino spectrum av- 
eraged over zenith angle [29] (see also [30]). 

Figure [4] presents the comparison of the calculation 
of the conventional (from fi, ir, i-T-decays) and prompt 
muon neutrino flux [H LH 131 131 El with the 
data of the AMANDA-II experiment g]. The conven- 
tional flux here was computed with use of QGSJET-II 
model combined with ZS primary spectrum (thin lines 
"conv."). Dashed lines mark the calculation by Nau- 
mov, Sinegovskaya and Sinegovsky [131 Q15] of the con- 
ventional muon neutrino fluxes for 6 = 0° and 90°. 
Bold dotted line (1) shows the sum of the prompt 




Figure 3: Conventional z/ M + flux at different zenith 
angles. Blue points: IceCube preliminary muon neutrino 
spectrum averaged over zenith angles [25] . 
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Figure 4: Fluxes of the conventional and prompt muon 
neutrinos. Experiments: AMANDA-II 4] (circles) and 
Frejus |3T] (the square). Calculations: the conventional 
flux - thin red lines (this work) and dashed [13]; the 
prompt flux - VZ [33] (line 1), RQPM [32] (2), GGV 
[34] (3, 5 for A = 0.5,0.1), QGSM [32] (4). 



neutrino flux by Volkova & Zatsepin (VZ) [33] and 
the conventional one due to the QGSJET-II + ZS 
model at 9 = 90°. Dash-dotted line (2) marks the 
sum of the QGSJET-II conventional flux (0 = 90°) 
and the prompt neutrino contribution due to the re- 
combination quark-parton model (RQPM) [31] . Solid 
line 4 shows the same for the prompt neutrino flux 
due to the quark-gluon string model (QGSM) [32] (see 
also [331 [19j [35] ). Also shown here are the two of the 
prompt neutrino flux predictions by Gelmini, Gondolo 
and Varieschi (GGV) [34] : line 3 (5) represents the 
case of A = 0.5 (0.1), where A is exponent of the gluon 
distribution at low Bjorken x. Curves just below 3, 
4, 5 ones display the coresponding flux at 9 = 0°. 
Calculated prompt neutrino fluxes at E v = 100 TeV 



are presented in Table III along with the upper limit 
on the astrophysical muon neutrino diffuse flux ob- 
tained in AMANDA-II experiment g]. Note that the 
QGSJET-II+GH flux appears to be the lowest flux of 
the conventional neutrinos at high energies. 



Table III Atmospheric neutrino flux at E v = 100 TeV vs. 
the AMANDA-II restriction for the + flux 



Model 


El4> u , (cn^ssr)" 1 GeV 


conventional + : 
QGSJET-II + ZS 
QGSJET-II + GH 


0° 90° 
1.20 x 10~ 8 10.5 x 10" 8 
1.11 x 10~ 8 9.89 x 10~ 8 


prompt + : 
QGSM [32] 
RQPM [3J 
VZ [33] 


90° 
1.22 x 10~ 8 
4.61 x 10~ 8 
8.12 x 10" 8 


AMANDA-II upper limit g] 


7.4 x 10~ 8 
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4. Summary 

The calculations of the high-energy atmospheric 
muon neutrino flux demonstrate rather weak depen- 
dence on the primary specrtum models in the energy 
range 10 — 10 5 GeV. However the picture appears 
less steady because of sizable flux differences origi- 
nated from the models of high-energy hadronic inter- 
actions. As it can be seen by the example of the mod- 
els QGSJET-II and SIBYLL 2.1, the major factor of 
the discrepancy in the conventional neutrino flux is 
the kaon production in nucleon-nucleus collisions. 

A common hope that atmospheric muon fluxes 
might be reliable tool to promote the discrimination 
between the hadronic interaction models seems to be 
rather illusive as the key differences in the tt/K pro- 
duction impact variously on the neutrino flux and 
muon one. For the high-energy neutrino production 
at the atmosphere the kaon yield in nucleon-nucleus 
interactions is more strong factor in comparison with 
that for production of the atmospheric muons, despite 
on their common to neutrinos origin. 

Inasmuch as the atmospheric prompt neutrino flux 
weakly depends on the zenith angle (near 100 TeV), 
one may refer the AMANDA-II restriction just to the 
prompt neutrino flux model. Thus one may consider 
both RQPM and QGSM to be consistent with the 
AMANDA-II upper limit for diffuse neutrino flux. 
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